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Abstract 

An individual defect in a nematic liquid crystal moves not only in response to its interaction 
with other defects but also in response to an external field. We analyze the motion of a wedge 
disclination in the presence of an applied field of strength H. We neglect backflow and seek 
steadily travelling patterns. The stationary picture yields a semi-infinite wall of strength n, 
bounded by the defect line. We find that the disclination advances into the region containing 
the wall at velocity v{H), where v scales as H/\\ogH\ as long as the magnetic coherence 
length is greater than the core radius. When the external field is applied in the presence of a 
pair of disclinations, their dynamics is strongly infiuenced. We compute the expected relative 
velocity of the disclinations as a function of distance and field. The natural tendency for the 
disclinations to annihilate each other can be overcome by a sufficiently strong field suitably 
directed. 



1 Introduction 



Singularities in liquid crystals, or defects, have played a fundamental role not only in the develop- 
ment of the understanding of the physics of liquid crystals but also in the later development of the 
topological theory of defects in condensed matter. In nematic liquid crystals, point, line and wall 
defects can be found. Line defects were first classified by F.C. Frank, who noted that line defects 
came in classes with an integer or half-integer charge 1 . Later workers, using the topological 
theory of defects, realised that nematic liquid crystals sustain a topologically unique line defect 
miSl^- However, Frank's naive classification, which effectively supposes that the nematic order 
parameter is restricted to a plane, remains important in providing guidance and intuition for the 
physics of defects in nematic liquid crystals. 

The topological total charge of a system is conserved during its evolution. Under many circum- 
stances this remains true for Frank's definition of charge; this is a stronger condition. For instance, 
defects of opposite charges may annihilate and defects of higher topological charges may decay to 
a bunch of defects of smaller topological charges. Topological dipoles may even be nucleated from 
a smooth field . 

In this paper we focus on nematic disclinations, i.e. line defects. According to Frank's defini- 
tion, the topological charge of a line defect is defined as the number of turns the director performs 
along a closed path surrounding the defect. This number may be half-integer because of the head- 
and-tail symmetry of nematic liquid crystal molecules. When the final director is rotated by an 
angle tt with respect to the initial director, the physical state they describe does not exhibit any 
discontinuity. 

The physical manifestation of the topological theory of defects in nematic liquid crystals con- 
cerns escape into the third dimension. Nematic disclinations of integer charge are not topologically 
stable. A suitable continuous transformation, involving the third previously neglected dimension, 
dissolves the singularity and leaves behind a regular field 6). In a similar manner, all half-integral 
disclinations can be distorted into each other. But often free energy criteria impede this pro- 
cess, and the Frank classification remains useful. In particular, annihilation of opposite charge 
half-integral defects is favored, but of same-charge defects is hindered. 

Defect dynamics has been widely studied either in the classical Oseen-Zocher-Frank (OZF) 
theory [3111^, or within the extended de Gennes-Ericksen theory [HlEIl- However, there are 
fundamental dilemmas in treating defect motion. The OZF theory is not a dynamical theory, and 
defect cores are regions in which, strictly speaking, the OZF theory does not apply. The theory can 
be extended to deal with moving defects, but only by introducing a phenomenological dissipation 
function using director rate of change. This minimal extension is incorrect; the full extension 
involves hydrodynamic terms which depend on the local director ^Jj . But the local director is not 
defined in the defect region, and so this extension is also inadequate to describe defect motion. The 
alternative de Gennes-Ericksen approach (now usually called the Q-tensor theory) is in principle 
up to the task, but now the defect regions are no longer special; there is a danger of confusing 
calculation with computation, leading to answers of doubtful physical significance. 

Despite these problems, some theoretical progress has been made, in part because some authors 
have detected analogous topological structures in liquid crystals and cosmological models. A single 
disclination may move through an otherwise smooth field 12 , but then the problem is why the 
defect is moving in the first place. An implicit response to this question is to give more attention 
to the interaction between two or more defects. Now the defects move, slowly, to reduce their 
elastic energy. For example, attraction between two point defects of opposite charges has been 
studied both in planar ^J^JE] and cylindrical |lfi[ 1171111^ geometries. The attraction between 
a dipole of smectic disclinations has been studied in . 

More precise quantitative descriptions of the defect evolution must necessarily take into account 
backflow effects [201 IS 122 j *-e- the interactions between director rotation and macroscopic 
molecular motion. The first analytical attempt to introduce backflow effects was performed in 
|28j ■ where the macroscopic velocity field was coupled to the degree of orientation, though not 
to the director field. A series of recent numerical simulations |^ I25L I26L 177] have determined 
the main effects of backflow coupling. The dynamical director patterns turn out to be strongly 
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influenced in the final part of the annihilation process, but some effects may be noted even during 
all the defect evolution. In particular, the positive-charge disclination moves almost twice as faster 
than the negative-charge one |25) . 

In this paper we investigate the effect of applied external field on the defect dynamics. We 
consider a simple geometry: a single or a dipole of ±i disclinations, in planar symmetry. 
We compute the defect speed with the aid of the Leslie's dissipation balance JI]. It turns out 
that the external field exerts a profound effect on the defect dynamics and the defect interaction. 
By suitably adjusting the external field strength and direction it is possible to drive a single 
disclination through the sample, almost at will. More interestingly, the external field may not only 
accelerate the annihilation process. It can also stop the defects, or even reverse their velocities, so 
transforming an attractive into a repulsive interaction. Depending on the external field direction, 
a critical defect distance may arise. This characterizes the defect interaction. They annihilate 
each other if they come closer than that distance; otherwise, they repel. 

Throughout our presentation we will introduce and discuss some assumptions that simplify our 
analytical calculations. The 1-constant approximation in the elastic free energy, and a parabolic 
approximation in the magnetic energy linearize the free energy derivatives. We remark that the 
parabolic approximation should be abandoned if we were to generalize the present study to highly 
charged defects. We also neglect backflow. This approximation allows us to derive an analytical 
condition which determines the defect velocity, and in particular the critical distance that reverses 
the defect interaction. A numerical analysis would correct these computed values, even if the 
described phenomena will certainly remain. 

The plan of the paper is as follows. In Section 2 we analyse the motion of a single disclination 
in an external field. Section 3 describes the defect interaction, and how it can be reversed with 
the aid of the applied field. Section 4 briefly summarizes and discusses our results. 

2 Single defect motion 

We consider a disclination embedded in an external magnetic field. This field will favor the 
director orientation of one side of the defect with respect to the other side. This asymmetry is 
sufficient, as we shall see, for energy considerations alone to determine that the defect will move, 
and to determine the direction of its motion. Our task is to determine the magnitude of the 
velocity as a function of the field intensity. We adopt the 1-constant approximation and neglect 
backflow. This latter approximation implies that our estimated velocity will be certainly smaller 
than the actual velocity. In fact, backflow effects reduce the total dissipation, so allowing faster 
director dynamics. 

Let us consider the reference frame illustrated in Figure ^ It moves with the disclination, 
which sits at the origin O. The a;-axis is parallel to the external field. Given any point P in the 
plane, let be the polar angle between the radial direction {P — O) and the x-axis. Disclinations 
of half-integer charge do not benefit from escaping into the third dimension. Thus, we restrict our 
interest to situation in which the director n is confined to the plane, and let ^(P) be the angle 
the nematic director at P determines with the z-axis. 

In the absence of a field, the 1-constant approximation implies that the director angle varies 
linearly with the polar angle: $ = ^ In the presence of the field the same linear dependence 
holds, roughly speaking, on length scales smaller than the magnetic coherence length ^. So long 
as we confine ourselves to these length scales, the elastic energy overwhelms the external field. 
However, far from the defect, the magnetic energy forces the system toward its preferred value 
<& = 0. On the other hand, over a closed path around the defect, no matter how far that path 
may be from the defect, topology forces the angle $ to rotate through tt. This creates a dilemma, 
for the system must rotate through an angle tt, and yet remain at $ = 0. 

These two constraints are reconciled by restricting the region over which $ rotates. This region 
is a topologically irreducible wall between regions of space whose director is oppositely oriented. 
Associated with the wall is a well-defined surface free energy analogous to the surface tensions 
of fluid mechanics. This is a 7r-wall, since the director concentrates its 7r-rotation in it. In the 
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Figure 1: Geometric setting for the analysis of the single defect motion. The defect occupies 
the origin of the co-moving reference frame, and its moves towards the 7r-wall. The angle $ is 
determined by the director and the motion direction (orthogonal to the external field) . 



presence of the field, the previously isolated defect line has been transformed into the trailing edge 
of a wall defect. The disclination then moves into the wall in order to reduce the wall area and 
consequently the free energy of the system. 

2.1 Dissipation principle 

We determine both the director field and the value of the disclination velocity by imposing the 
energy balance between the free-energy loss-rate and the dissipation 

T + V^O. (2.1) 

In the 1-constant approximation, the free-energy functional is given by 

= / fif|V$P +Xai?^sin2$) da , ^2.2) 

where K is an average elastic constant, Xa is the magnetic anisotropy, and H is the strength of 
an external magnetic field. 

It is well-known that a -|-^ disclination possesses an infinite core energy. The elastic free-energy 
density diverges at the defect, and the singularity is not integrable. There are two ways to avoid 
this divergence. The first consists in excluding from our integration region a small disk centered 
in the defect, the core region Bo- The radius of the excised disk, the core radius Tq, is usually 
much smaller than all other characteristic lengths entering the problem, so that many studies 
have been performed in the limit of vanishing To |28j . A more precise physical description of the 
defect requires an extension of the classical Oseen-Zocher- Frank theory, and the replacement of the 
nematic director with the nematic order tensor [29) . We will choose the first option, and perform 
all the integrations in the pierced domain, which excludes Bo- We further assume that the core 
radius is fixed. The basic physics of the phenomenon is well-described using these approximations. 
However, in order to deal with external fields of any intensity, it would be interesting to apply the 
techniques of |3U| to determine how and when the magnetic coherent length influences the core 
radius. 

A second, though less worrying, free-energy divergence comes from the supposedly infinite 
size of the domain. The domain may extend indefinitely in the y-direction without inducing any 
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singularity, since both the elastic and the magnetic energy densities vanish away from the x-axis. 
On the contrary, there will be few cases in which a large-x cut will be needed to keep energies 
finite. In those cases we will assume that our domain is bounded by |x| < L. Whenever possible, 
we will perform the L — > +00 limit, and we will notice that the large scale length L will not finally 
enter in the defect velocity. 

Our final approximation concerns the magnetic free-energy. We will perform the parabolic 
approximation sin^ $ ~ This approximation allows us to obtain linear field equations in 
It can be used in the whole domain provided we define $ € [~f 1 f ); since it that case all the 
values attained by the director angle belong to the potential well of the equilibrium configuration 
$ = 0. We remark that this approximation would not be valid if we were interested in analysing 
the equilibrium configurations of more complex defects. 

When we neglect backflow, the dissipation assumes the simple expression 

V = ji ^^dxdy , (2.3) 
J Bo 

where the pierced integration domain Bo comes into play since also the dissipation density diverges 
in the core region. 

Let us perform an infinitesimal displacement (5$ of the director field. If we make use of the 
divergence theorem, the dissipation principle H2.1|l requires 

/ 5<^ {2KW<S>) ■ d£+ (5$ [71 $ - 2XA'I> + 2xa-ff^$j da = V^*, (2.4) 

where v is the outer normal along dBo. The arbitrariness of 5$ requires that the quantity in 
square brackets in the second integral must vanish identically: 

71 $ = 2XA$ - 2xai?^ * • (2.5) 

Equation H2.5|l is the well-known time-dependent Ginzburg-Landau evolution equation of the sys- 
tem. It will supply us the director field. However, the Ginzburg-Landau equation alone is not able 
to guarantee the dissipation principle. We must check that also the first integral in (|2.4|l vanishes. 
The boundary dBo is made of two parts: a small circle around the defect and a large boundary 
at infinity. The integral around the latter vanishes since at infinity V*!" fades everywhere, except 
along the 7r-wall, where it is orthogonal to the outward normal. We are then left with the first 
integral in 12.4|l . performed along the boundary of the core region. This quantity has a simple 
physical meaning |31j : it is the power supplied by the core region to the rest of the domain. Thus, 
to require that all the free-energy loss is dissipated within the system is equivalent to require that 
no energy is supplied to it, neither from the outside nor from the core region. 

2.2 Steadily moving defects 

We will look for stationary solutions of equation H2.5(l . They aim at representing a defect moving 
at a constant speed v towards the 7r-wall. We will find that, for any positive value of v, it is 
possible to find a solution of 1)2. 5|l which satisfies the boundary conditions. However, there is just 
one value of the velocity which cancels also the first integral in l|2.4(l (or, equivalently, that satisfies 
the global dissipation condition). It is the only velocity at which the defect is able to move without 
any external boost. 

In a steadily moving reference frame, traveling with the defect itself, the relative stationary 
differential equation follows by replacing in (|2.5|l the time derivative with vdx- We thus obtain: 

A 9$ $ 



Ik. . V , , , 

where £, ■= \ tt is the magnetic coherence length, and A := — ^= - — r is a dimension- 

V Xai?^ 2y/KXe. \H\ 

less quantity. We look for solutions of the eigenvalue problem ()2.6|l that satisfy the symmetry 
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requirement ^{x, —y) = 



— y) for all y ^ 0, and the boundary conditions 



$(a;,0+) 



if x > 

- if X < , 
2 



and 



lim '^(x, y) = for all x e 

J/— »oo 



(2.7) 



where the latter condition is determined by the presence of the magnetic field. The boundary 
conditions H2.7|l are singular only at the origin, where a disclination of topological charge +5 
stands. Indeed, the discontinuity the angle <f> suffers along the negative x-axis does not induce 
any physical singularity, since $ = 7r/2 and $ = — 7r/2 describe the same director orientation. 

Among the solutions of the eigenvalue problem (|2.6|) - (|2.7|l . the dissipation principle (|2.1() will 
single out the only physical one. 



2.3 Director field 

We solve the linear partial differential equation H2.6|l with a Fourier transform. To this aim, it is 
useful to write the first of the boundary conditions H2.7|l as 



$(a;,0^ 



TT 1 

4-41 



dq 

q 



for all X ^ 



(2.8) 



where PV denotes the Cauchy Principal Value of an integral. We then look for solutions of 
equation H2.6|l of the form: 



with 



^{x,y) = -9i{y) - - FYj^^g2{q,y)dq 



ffi(O) = 92{q, 0) = 1 and lim gi{y) = lim 52(9, y) = for aU q e 



If we insert (jTU in we obtain: 



9'{ 



■91 



41 



PV 



q 



d'92 

9y2 



kHq)92 



dq^Q, 



where k{q) will henceforth denote the positive-real-part square-root of fc^(g) = 
The solution of H2.11|l and 1)2.10(1 in the upper half-plane {y > 0} is 

l^iqx-k{q)y 



iqX 

T 



■n _y \ f c^i^^'^wy 

^{x,y) = - e « - - PV/ dq , 

4 4i 7r g 



whereas, by symmetry, the solution in the whole plane is given by: 



<^{x,y) = e{y) 



1 

4i 



PV 



^\qx-k{q)\y\ 



■ dq 



(2.9) 
(2.10) 

(2.11) 
1 

(2.12) 
(2.13) 



with 



(2.14) 



1 ify>0, 
^-1 ify<0. 

Figure [21 illustrates the solution (|2.13|) . More precisely, in the upper plot it is possible to observe 
how the magnetic free-energy density is concentrated on the 7r-wall, whereas the lower plot shows 
that the elastic free-energy density is mostly localized on the defect. In both cases, it is possible 
to notice that the decay pattern from $ = ±-| on the left a;- axis to the equilibrium value $ = 
becomes constant when we move some magnetic coherence lengths away from the defect. 

We notice that no use of the core region has been made for the time being. In fact, the 
differential problem l|2.6|l - H2.7|l admits a solution in \ (0,0), without any need of excising a 
finite region around the defect. We will see below that this will not be the case when we will have 
to deal with the derivatives of the field, and in particular with the dissipation H2.3|l . 
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Figure 2: Density plots of the magnetic (upper graph) and elastic (lower) parts of the free-energy 
density, computed from the analytical solution H2.13|l . The elastic energy is localized close to the 
defect. On the contrary, the magnetic energy is mostly packed around the 7r-wall. The defect core 
is not symmetric with respect to the defect position. It is slightly shifted behind it, with respect 
to the direction of motion. 



2.4 Disclination velocity 

We have succeeded in finding a solution of the stationary Ginzburg-Landau equation for any value 
of V. We will now complete the eigenvalue analysis, and determine the correct disclination velocity, 
by imposing the global dissipation condition H2.1|l . 

We have already noticed that, in the moving reference frame, time derivatives transform into 
spatial derivatives. Thus, 



d_ 

dx 



X — + 00 



(2.15) 



When |a;| is very large, the principal value of the g-integral in H2.13|l is dominated by the low q 
values: 



lim X PV 

:^±oo 4i 



■dq = X lim PV 

4i X — »-±CXD 



Thus, 



7" = 



e 



tt'^Kv /■+°° 



e ^ dy — 



q 

n^Kv 



dq = ±^ e-«/« 



(2.16) 
(2.17) 
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The result H2.17|l admits a simple physical interpretation, that aheady exhibits in H2.15|l . The 
quantity within square brackets in (|2.15|l is the free-energy contained in an infinite vertical strip 
of unit width, centered at x. When x — > +00, both the elastic and magnetic energy densities 
relax to 0, as Figure |21 shows. As a consequence, the total energy stored in the right-side strip 
is negligible. The picture completely changes in the x —00 limit. A unit-width strip crossing 
the TT-wall stores a finite amount of free energy, measured in H2.17(l . In fact, it is precisely the 
difference between the energies stored in those strips that keeps the defect moving. In the unit 
time, the defect motion replaces a left-side strip of width v, with free energy given in (|2.17|l . with 
a right-side strip of equal width, with no free energy. 

While the free energy variation depends only on the asymptotic structure of the director field, 
dissipation takes place in the whole domain. Indeed, it is so strong close to the moving defect, 
that we will be forced to exclude the core region in order to avoid an infinite dissipation. We have: 



X> = 71 / $ (ix = 71 u / ( -g— ) dx dy 



^ f dxdy Pv/'dg Pv/'d<z'e'(«+^')-('=(«)+'=(9'))l^'l 
PV/ dq PV/ dq' ^ ' 



k{q) + k{q') 
7r7iw2 /■'?M/(27r) 



2 7o v?+We+W+vF^^We+W 



2^ Jo 




(2.18) 



where 5 denotes the Dirac delta function. The high-q cutoff is needed in order to avoid the 
logarithmic divergence which the disclination induces both in the free energy (but not in its time- 
derivative) and in the dissipation. This is related to the inverse of the core radius: qyi = 27r/ro. 

When we replace (|2.15(l and H2.18|l into the dissipation principle H2.1|l we obtain the self- 
consistency equation that determines A (i.e., f), as a function of the ratio i,/ro'. 

(2.19) 

The most interesting region in physical applications is ^ 3> ^o- The integral on the left hand side 
is dominated by its logarithmic high-s divergence and we obtain: 




iAlog — = i^ =^ v= \H\ when ^ > • (2.20) 

To 71 log (^/r-o ) 

However, it is interesting to push the analysis of H2.19|l into the opposite regime ^ <C Tq. The low-s 
terms dominate the integral on the left hand side and we find: 

iA— ^iTT =^ v^^^^MIh^ when ^ < f^o . (2.21) 

fo 71 

The quite intriguing quadratic behaviour predicted by 12.21|l must be handled carefully. When 
the external field becomes so intense that ^ becomes of the order of ro, we have to question 
our assumption of ro being independent from ^. A more complete theory, that can be derived 
following the steps of |30| . would yield an ro(^), and thus a disclination velocity depending only 
on the strength of the applied field. 

Figurcinishows the numerical solution of equation (|2.19|) . The plot highlights that the transition 
between the two asymptotic regimes derived above occurs for ^ > ro. Even though this is the 
limit up to which we can seriously trust our analytical result. Figure |31 suggests that the nonlinear 
effects increase the disclination velocity. 
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Figure 3: Disclination velocity as a function of the external applied field 

3 Defect interaction 

Let us now consider two nematic defects of opposite topological charges s = ±i placed on the 
X-axis at a distance A — 2d, but in the absence of any external field. After a certain period of 
time these defects amalgamate and thus annihilate each other. The two cores disappear at once, 
with a large and rapid reduction in free-energy |13l 1151 1^ . Imagine, however, that we apply a 
magnetic field that favours the director orientation of the molecules that lie within the defects. 
The defect speed will be certainly reduced. In this section we show that it can be even reversed. It 
exists a critical distance such that the defect interaction is attractive only if their mutual distance 
is smaller than the critical one. Otherwise, they repel. 

Throughout this section, we will work out in detail the geometry illustrated in Figure 01 The 
applied field lies parallel both to the line connecting the defects, and to the director orientation 
between them. In this geometry, the defects may only approach or separate, so avoiding but more 
complicated motions such as mutual rotations. 

We compute the speed of two stationarily moving defects. We will show below that the sta- 
tionarity assumption holds approximately even when the defect distance becomes of the order of, 
or smaller than, the magnetic coherence length. It must certainly be abandoned when one wants 
to describe the complete annihilation process, and in particular when the defect distance becomes 
of the order of the core radius. 

In the final subsection we will briefly analyse the case when the external field determines 
a generic angle a with respect to the defect line. In this case the motion may be much more 
complicated. However, we are able to estimate how the nature of the defect interaction depends 
on a. More precisely, we will determine for which values of a the defect interaction may become 
repulsive, and how the critical distance depends on the external field direction. 

The time-dependent Ginzburg-Landau equation H2.6() is linear, due to the parabolic approxi- 
mation we used for the magnetic energy. We can thus obtain a solution describing two stationarily 
approaching (or separating) defects by simply superposing two functions of the type (|2.13|l . More 
precisely, we add a solution describing a defect placed at x = +d, travelling with velocity —v, and 
a solution describing a defect placed at x = —d, travelling with velocity The velocity v{d) 
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Figure 4: Two attracting defects placed at a distance A along the x-axis. The external field is 
parallel both to the line connecting the defects and to the director orientation between them. 



will be determined later, again using the dissipation principle. The director field is given by: 
^ix,y) = e{y) -e ^ 



£(2/) 



1 /■ pi9(2^-<i)-fc('?)|y| _ pi«(3:+d)-fc(-g)|i/| 

1 PV/ dq 

4i Jr q 



■e « 



i qx sin {qd + Im k{q) \y\) e" 



dq 



(3.1) 



Figure El illustrates the countervailing tendencies of the elastic and magnetic energies: the 
elastic contribution aims at annihilating the defects in order to relax the infinite core energy. On 
the contrary, the magnetic field tries to broaden the intermediate region, where all the molecules 
are already correctly aligned. 

The free energy associated with the traveling configuration (I3.1|l is given by: 



J^(A, A) = Ktt 



h arcsmh A — 



ttA 



sm 



2^ 



ds + Fi{A, X) + F2{X) 



(3.2) 



The first term in l|3.2() diverges in the L +00 limit {L is the horizontal scale of the system: 
X € [—L, L]). This term corresponds to the energy of the two walls. The second term is related to 
the infinite elastic energy stored in the defects (ro is the core radius as in the preceding section). 
The next three terms depend only on the defect distance, but not on their velocity. In particular, 
the term containing the modified Bessel function of the second kind Kq is able to cancel the free 
energy divergence when A ~ Tq. Finally, both Fi and F2 are velocity corrections, which vanish if 
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Figure 5: Density plot of the magnetic free-energy density (upper graph) and the elastic free- 
energy density (lower graph), computed from the analytical solution H3.1|l . For these particular 
values of the distance and the magnetic coherence length, it will turn out that the defects are 
separating (that is, v < 0). Both the magnetic and the elastic free-energy densities are mainly 
concentrated on the defect walls. Between the defects, the free-energy density is negligible. 



A (i.e., v) vanishes: 
F,iA,X) 

F2{X) := 



1 

2 7o 



Re 
Re 



^I,{\,s) - m(o, s) + 
1 1 



/i(A,s) m(0,s) 



1 1 



ds 



1 



/i(A,.s)+/i(-A,.s) 2/i(0,s) 



ds 



1 



-isA/C 



T2 



ds 

72 



71 



(3.3) 



(3.4) 



4^0 m(A, s) -I- /x(— A, s) \ fi{X, s) fi{— X, s) 
where /i(s) is the positive- real-part square- root of /i^(s) = 1 + is A + s^, and A „„ 

SV-ft^Xa \H\ 

in the preceding section. 

In the absence of backflow, the dissipation stems only from the director rotation. When 
computing $, we assume that only the defect position d depends on time. In our stationary 
approximation, we thus neglect the time derivative of the velocity v. However, we remark that 
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when we want to determine the non-stationary effects depending on i), we are no longer aUowed 
to use the solution (|3.1|) . The dissipation function is given by: 



I? = 71 



with 



dx dy = 
k2 2 



Re 



Gi(A,A) := 

G2(A) := 2 / 

"'0 



arcsinh ■ 



1 



1 



+ Gi(A,A) + G2(A) 



^(A, s) Ai(0, s) 



Ai(A,s) +^(-A,s) 2//(0, s) 



ds 
ds 



and 



(3.5) 

(3.6) 
(3.7) 



The functions Gi and G2 vanish in the low-velocity limit v 0. The dissipation principle delivers 
the self-consistency equation that determines the defect velocity v = —d: 



'dd 



d + V = 



A 



(3.8) 



^ + ^ 



sin s ds 



A2 



,2 ^ 9A 2 



A 



arcsinh + /sTq ( y ) + Gi ( A, A) -|- G2(A) 



Before analyzing the solutions A(A) of equation H3.8|l . we want to stress the importance of the 
first two terms in its left-hand side. In fact, if we define 



fix) ■.^Ki{x)+x 



sm s 



we have: 



5^ 
OA 



v=0 



Ktt 
T 



sV x^ + 



(f 



ds 



(3.9) 



(3.10) 



Thus, in general, the defects will approach or separate depending on whether /(A/f) exceeds tt 
or not. We postpone the analysis of the properties of / to the next subsection, when we will 
generalize equation (|3.1()|l to the case of tilted applied fields. As for now, we only remark that 
f{x) = TT when x = Xcr = 0.377388. The function / is greater than tt (thus inducing defect 
attraction) when A < XcrS.- Defect repulsion is induced at distances greater than cCcrC- 

Figure illustrates the numerical solutions of equation H3.8|) for three different values of the 
ratio between the magnetic coherence length and the core radius. They exhibit the following 
properties. 

• In the large distance limit, all Bessel functions decay exponentially with A/^. The functions 
Gi, G2, and the derivative of Fi vanish too. Furthermore, 



lim X 

X — >oo 



sm s 







s^fxP' + s^ 

Thus, the large-distance limit of A is given by 



ds 



sm S TT 

ds = — 

s 2 



lim A(A) 

A^oo 



arcsinh(^/ro) 



V 



27rV7?^ \H\ 
71 arcsinh (^/ro) 



(3.11) 



whenA>C. (3.12) 



The defects repel, and move at a constant speed. In fact, if we compare H3.12|l with (|2.20|l . 
we notice that in the large-distance limit the defects behave independently, each moving at 
the velocity computed in the 1-defect case, since arcsinh(^/ro) ~ log(^/ro) when ^ ^ Tq. 

• The critical distance at which the defect interaction changes sign does not depend on the 
core radius Vo'- all plots in Figure El cross the x-axis at A = Xcr^- 
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• In the small distance limit, the stationary approximation we have used is not well justified. 
In this limit, the velocity diverges, and all terms depending on v must be taken into account 
in the equation of motion. In any case, we remark that the Bessel function Ki induces a 
divergence in the derivative of the free-energy that scales as the inverse of the defect distance. 
Figure El suggests that the fully nonlinear regime is limited to distances much smaller than 
the magnetic coherence length. 

10 IT 1 




Figure 6: Defect velocity as a function of the external field for different values of the core radius 
ro- C/^o = 10 (dotted fine), 100 (dashed), 1000 (full). For any value of the core radius, the defects 
attract only if their distance is smaller than Ad-. 



3.1 Tilted external fields 



Let us now imagine that the external field is rotated of an angle a with respect to the direction 
of Figure ^ In this case, the nematic director is forced to relax to the external field direction 
everywhere but on the topologically irreducible 7r-wall. The motion will be much more complex 
- the defects will rotate, and the 7r-wall will not necessarily be straight at all times. However, 
it is possible to ascertain whether the defect interaction will lead to attraction or repulsion. To 
this aim, we imagine to pin the defects at a distance A, and we look for the stationary director 
field in the presence of a tilted external field. Then, we compute the free-energy of the stationary 
solution, and we check the sign of its derivative with respect to the distance. We stress that it is 
not possible to use this derivative in a dissipation principle to obtain a defect velocity. However, 
its sign will determine whether or not the defects, whatever their complex motion, will approach 
each other. 

The stationary director field can be simply derived by running through the above steps again. 
In the presence of two defects placed at a; = ±d, and an external field tilted of an angle a with 
respect to the a;-axis, the stationary configuration is given by: 



*(a;,y) = £(y) 







Ivl f 


a + 




e-e - / 






Jo 



cos qx sin qd e 



-ko(q) \y\ 



dq 

q 



(3.13) 



where ko{q) coincides with the zero-velocity limit of k{q) above: ko{q) :— \/ q'^ + 1/^^ ■ 

If we compute the free-energy associated to H3.13|l . and we differentiate it with respect to A, 
we obtain the generalization of equation (|3.10(l : 



dA 



Ktt 










;!)-(.- 2.); 


■u=0 ^ 





with the same / defined in 
following properties: 



(3.14) 



The left panel of Figure [T] shows the plot of /. It enjoys the 
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• f'{x) = i {^Ko{x) — K2{x)^ < for all x > 0. Thus, the critical distance at which the defect 
interaction becomes repulsive increases when the external field is tilted. Furthermore, the 
free energy is a concave function of the distance between the defects. 

• lim f{x) = +00 and lim f{x) = — . More precisely, 

f{x) = - - ^ logx + 0{x) as a; ^ 0+ and f{x) = '^+0 (s"^) as a; ^ +00 . (3.15) 
Thus, the equation 

dT 

— = ^ /(a;)=7r-2a, (3.16) 

which determines the equilibrium distances, possesses exactly one solution if and only if 
a < J . This limiting value for the tilting angle could be easily predicted. In fact, if 
the external field determines an angle greater than j with the defect line, the director 
orientation between the defects costs more magnetic energy than the director orientation 
outside the defects. Thus, in this case, the external field strengthens the defect attraction 
at any distance. 

In summary: 

• When a = 0, the defects attract if A < Acr(f) = 0.3774^; they repel when A > Acr(f ). 

• When a > 0, the critical distance Acr(/3) increases. It diverges when a — > 

• When a > J there is no critical distance: the defects always attract (the external field 
enhances the attraction). 




Figure 7: Left: Plot of the function / defined in the text. Right: Critical value of the defect 
distance as a function of the tilt angle a. 



4 Discussion 

We have studied the motion of a single disclination, and a disclination dipole in an external field. 
Our results show that it is possible to drive the disclinations by suitably adjusting the external 
field direction and strcmgth. 

In the case of a single disclination we have shown that the disclination velocity depends almost 
linearly on the field strength, since the coefficient depends on the logarithm of the field strength 



13 



(see equation (|2.20|l ). Figure |21 shows that the hnear scaling is abandoned when ^ ~ Tq. Equation 
(|2.21|) shows that below this regime the disclination velocity is expected to scale as the square of 
the applied strength. However, this last prediction should be tested carefully, since it is surely 
influenced by our assumption that the core radius does not depend on the field strength. 

We have also computed the velocity of two opposite disclinations in the presence of a field that 
promotes the director orientation between them. The velocity depends on the defect distance A 
but, whatever the value of the ratio ^/ro, they defects approach if A is smaller that a critical 
distance A„] otherwise, they repel. We have finally generalized our calculations in order to deal 
with rotated external fields. Figure [T] (right) shows how the Acr depends on the angle a that the 
external field determines with the lie which connects the defects. 

The introduction of backflow may change the picture we have developed, both in the one and 
two-disclination cases. In the one disclination case the dissipation may sometimes be significantly 
reduced. The system can match the disclination motion with a flow that almost cancels the 
dissipation in the crucial core region. Whether this occurs seems likely to depend on the charge 
of the central disclination. We note that central physical and mathematical issues associated with 
backflow in nematic liquid crystal problems still remain open even over a quarter of a century 
after its essential role was first realized '20'. 

Likewise in the two-disclination problem, our computed defect velocity is symmetrical in both 
disclinations. This is also an effect of the no-backflow simplification. A generalization of the 
present study should allow to compute the different approaching or separating velocities. It could 
even happen that one defect moves towards the other, but that the other retreats faster still, 
allowing the defects eventually to separate. However, the existence of a critical distance that 
reverses the defect interaction can not be erased by backflow effects. The critical distance stems 
from the balance between the elastic and magnetic free-energy gains. The magnetic gain does not 
depend on the defect distance, while the elastic gain vanishes when the defects move apart. Thus, 
at some intermediate distance one overwhelms the other. 
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